Coherent Phonon Coupling to Individual Bloch States in Photoexcited Bismuth 
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We investigate the temporal evolution of the electronic states at the bismuth (111) surface by 
means of time and angle resolved photoelectron spectroscopy. The binding energy of bulk-like 
bands oscillates with the frequency of the A\ g phonon mode whereas surface states are insensitive 
to the coherent displacement of the lattice. A strong dependence of the oscillation amplitude on the 
electronic wavevector is correctly reproduced by ab initio calculations of electron-phonon coupling. 
Besides these oscillations, all the electronic states also display a photoinduced shift towards higher 
binding energy whose dynamics follows the evolution of the electronic temperature. 



PACS numbers: 

It is well known that most elements crystallize in struc- 
tures with hexagonal or cubic symmetry. One notable ex- 
ception is bismuth, which instead crystallizes in the A7 
rhombohedral structure. It has been originally proposed 
by Jones and Peierls that cubic bismuth would be unsta- 
ble towards a lattice distortion leading to the doubling 
of the unit cell pQ. Due to the distorted nature of the 
groundstate, the atomic positions are very sensitive to 
the electronic distribution in the conduction band. This 
conjecture has been experimentally verified by suddenly 
changing the occupation number of the electrons via the 
absorption of a femtosecond laser pulse [2H3] . Due to 
photoexcitation, a large coherent phonon modulates the 
distance between the two atoms in the unit cell along 
the (111) direction. The displacement vector indicates 
that the lattice coordinates move towards a hypothet- 
ical phase with a single atom per unit cell. In reality 
the material displays a non-thermal melting before the 
symmetric phase can be attained |5|. 

The atomic motion following the photoexcitation of 
bismuth is nowadays well understood. Time resolved X- 
ray diffraction and Density Functional Theory (DFT) cal- 
culations could accurately describe the amplitude of the 
oscillations and the frequency softening at high excita- 
tion densities [BHB]. On the other hand, the temporal 
evolution of electronic states has never been directly ob- 
served. Some aspects of the electronic dynamics can be 
inferred from general arguments. The generation mech- 
anism of coherent phonons with large amplitude implies 
that low energy electronic excitations are coupled to the 
Ai g mode. It is also expected that the electron-phonon 
coupling of such distorted material varies with the elec- 
tronic wavevector. Nonetheless, the determination of the 
coupling matrix elements for individual Bloch states is 
still an experimental challenge. We approach this task by 
time resolved photoelectron spectroscopy measurements. 
The electronic states of the bismuth (111) surface are 
characterized by angle resolved photoemission and first 
principles calculations. We confirm that the electronic 



structure displays a rich combination of bulk-like bands, 
surface states and surface resonances [9j[TT]. Upon pho- 
toexcitation, these electronic states display a dynamics 
that depends on wavevector and band index. The con- 
nection of our data with DFT calculations provides a 
full characterization of the coupling between electrons 
and the A\ g mode. In addition, we observe that the 
electronic states do not follow rigidly the motion of the 
lattice displacement. Instead, the bands also display a 
purely electronic shift towards higher binding energy. 

The experiments were performed with the Fem- 
toARPES setup, using a Ti: Sapphire laser that generates 
35 fs pulses centered at 810 nm with repetition rate of 
250 kHz. Part of the beam is employed to generate the 
fourth harmonic by frequency doubling in BBO crystals 
(/3-BaB 2 04). The 205 nm probe and the 810 nm pump 
are focused on the sample with a spot diameter of 100 
/xm and 200 /im, respectively. Their cross-correlation in a 
BBO crystal has a full width at half maximum (FWHM) 
of 80 fs. An electrostatic spectrometer with energy reso- 
lution better than 10 meV and angular resolution better 
than 0.5 degrees analyses the photoelectrons emitted by 
the 205 nm beam. The (111) surface of bismuth has been 
obtained by sputtering and annealing cycles of a single 
crystal. All measurements have been performed at the 
base temperature of 130 K and at the base pressure of 
7-10 -11 mbar. The out-of-equilibrium spectra have been 
collected with incident pump fluence of 0.6 mJ/cm 2 . 

Figure [T]A. shows an intensity map of the photoelec- 
trons emitted along the T-M direction of the Brillouin 
zone. As previously reported, the break-down of transla- 
tional symmetry in the (111) direction generates surface 
states that intersect the Fermi level and give rise to a 
Fermi surface [9 . These evanescent wavefunctions are 
localized at the topmost layers, thus conferring to the 
surface of bismuth good metallic properties. Although, 
the band structure supports surface states only near the 
Fermi level, some bulk-like states and surface resonances 
are clearly visible at higher binding energy [10] . Two 
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FIG. 1: A): Intensity map of photoelectrons emitted for par- 
allel wavevector along the T-M direction. Band a and band 
b are the surface resonance and the bulk-like state, respec- 
tively. The blue line is a bulk band calculated for perpen- 
dicular wavevector k± = 0.45 A -1 . B): Intensity map of 
photoelectrons emitted at binding energy e = —0.26 eV. An 
intense structure originating from band b is visible along the 
T-M direction but not along T-M'. 



bands of different nature are indicated in Fig. [T]A by let- 
ter a and b. Due to the matrix elements of the photoemis- 
sion process, band b is clearly visible only when the par- 
allel component of the electronic wavevector is ku < 0.25 
A -1 . Band a has larger binding energy, disperses more 
strongly and intersects band b at the T point. The overall 
agreement between our measurements and one step pho- 
toemission calculations is remarkable [11]. A. Kimura et 
al. [11 also show that band a has a large spin polar- 
ization whereas band b does not. This finding suggests 
that a is a surface resonance whereas b is a bulk band. 
In order to confirm it, we performed DFT calculations of 
bulk bismuth within the framework of the local density 
approximation and the ABINIT [12] package. Since the 
photoemission process does not conserve the perpendic- 
ular component of the electronic wavevector, the connec- 
tion between a peak in the photocurrent and a calculated 
band is formally not allowed [13]. Nonetheless, Fig. [T|\ 
shows that the band calculated along the T-M direction 
for perpendicular wavevector k± = 0.45 A -1 matches 
reasonably well the experimental intensity distribution. 
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FIG. 2: A): Photoelectron spectra acquired 100 fs before and 
300 fs after the arrival of the pump beam along the T-M di- 
rection at k\\ = 0.12 A -1 . B): The effective temperature of 
photoexcited electrons as a function of pump probe delay ex- 
tracted from the data (marks) and simulated with the model 
of ref. [8] (red line). C): Intensity map of photoelectrons 
emitted at wavevector k\\ =0.12 A -1 along the T-M direction 
as a function of pump probe delay. 



More insights on the dispersion of band b can be obtained 
by mapping in the reciprocal space the photoelectrons at 
fixed kinetic energy. Figure shows an intensity map 
of photoelectrons acquired at different wavevectors for 
binding energy e = —260 meV. The structure b is very 
intense along T-M whereas it vanishes in the direction 
T-M'. The experimental evidence that band b is invari- 
ant for a rotation of 120° but not for a rotation of 60° 
corroborates the bulk character of band b [9 . 

Photoexcitation by an intense and ultrafast laser pulse 
generates sizable effects on the photoelectron current (see 
movie in the supplementary material). Figure [2]A shows 
the spectrum acquired at fcy = 0.12 A -1 along the T- 
M direction for two different pump-probe delays. Note 
the large reduction of spectral weight taking place 300 fs 
after the arrival of the pump beam. The partial deple- 
tion of states a and b can be approximately described by 
a Fermi-Dirac distribution with an effective temperature 
of 1900 K. As shown by Fig. [2j3, the electronic temper- 
ature decays as a function of pump probe delay with a 
time constant of 6 ps. The relaxation of the electronic 
energy density at the surface is due to the combined ef- 
fect of electronic diffusion and heat dissipation into the 
lattice bath. We could accurately reproduce the tem- 
poral evolution of the electronic temperature by solving 
the coupled differential equations in ref. [8] with incident 
pump fluence of 0.6 mJ/cm 2 , electron thermal conductiv- 
ity ko = 0.04 W m _1 K _1 and a temperature dependent 
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electron-phonon coupling. 

Beside the large changes of the electronic distribution, 
we also observe strong variations in the structure of the 
electronic states. As shown by Figure [2]A., both peak a 
and b shift towards higher binding energies after the ar- 
rival of the pump pulse. In addition, the position of peak 
b displays periodic modulations at the frequency of the 
A\ g mode. Such a dynamics arises from the coupling 
of the coherent phonon to electronic states and can be 
directly observed on the photoelectron intensity map of 
Fig. [2p. The oscillations are large for the bulk band b 
whereas they fall below the detection limit for the sur- 
face resonance a and for all other surface states. This 
observation leads to two possible conclusions: either the 
amplitude of the A\ g mode is smaller on the topmost bi- 
layer or the surface states have a weak coupling to this 
mode. 

The excitation of coherent phonons in semimetals has 
been the subject of many experimental and theoretical 
works [2HH] • Since bismuth absorbs photons at 810 nm, 
the generation mechanism is mainly due to a displacive 
excitation of the charge density [2j |3]. The DFT cal- 
culations can help us to illustrate this point: figure [3]A. 
displays a colorscale plot of the electronic density p(r, Ti) 
at Ti — 130 K in the (110) plane. The sudden increase 
of electronic temperature from 130 K to Tf = 2000 K 
changes the spatial ditribution of electronic density. We 
show in Fig. [3^ that (p(r,l» - p(r,7-))/ p(r,Ti) is neg- 
ative in the region between nearest neighbours (atoms 
1 and 2). The transfer of electronic charge out of the 
bonding region increases the repulsion between the ion 
cores. As a consequence, the atoms oscillate around a 
new equilibrium position until the energy dissipation into 
other phonon modes brings the system back to the initial 
state. The lattice motion changes the binding energy of 
the Bloch state |6, fey) by A?,/^ ■ u, where u is the Ai g 
phonon coordinate and is the deformation poten- 

tial. Figure [3p shows the binding energy of peak b at 
the base temperature of 130 K as a function of pump 
probe delay. The oscillations induced by the A\ g mode 
have frequency v — 2.97 ± 0.05 THz and damping time 
7 = 2.6 ± 0.2 ps. These values compares well with the 
frequency and damping time observed by transient re- 
flectivity in comparable experimental conditions [4J. We 
reproduce the oscillations of peak b by solving the model 
proposed in Ref. [8 and using the value of the computed 
deformation potential A?,^ = 0.91 eV/A for fey = 0.12 
A -1 . The simulated lattice dynamics attains a maximum 
displacement value u = 1.6 pm, corresponding to 0.14 % 
of the body diagonal length c = 11.86 A. 

As shown in Fig. |3p, the measured amplitude of oscil- 
lation is 20 meV for fcy < 0.05 A -1 and decreases down 
to 10 meV for fcy < 0.2 A -1 . At larger wavevectors the 
photoelectron intensity of band b is too low for a reliable 
extraction of the D^ku ' u value. The large fcy -dependence 
of A>,/c|| has been previously reported for the amplitude 
mode of TbTe3 [15]. Due to the presence of a Charge 
Density Wave, the variation of the interaction strength 
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FIG. 3: A): Colorscale plot of the electronic density at 130 
K in the (110) plane. B): Relative change of the electronic 
density for an increase of electronic temperature from 130 K 
to 2000 K. C): Binding energy of peak b for k\\ = 0.12 A -1 
along the T-M direction as a function of pump probe delay 
(blue marks) and temporal evolution obtained by ab initio 
calculations (red line) [5]. D): Oscillation amplitude of peak 
b as a function of h\\ (marks) together with the prediction 
of DFT theory for a lattice displacement u = 1.6 pm (solid 
line). E): Calculated energy dispersion of band b along the 
T-M direction when the A\ g phonon mode is activated (u ^ 
0). The calculations are performed for two atoms located 
at ±(0.2329 • c + u) along the diagonal of the unit cell for 
u/c = 0.1 %, 0.2 % and 0.3 %. 



as a function of the bloch state has been explained by 
the good nesting properties of the Fermi surface. On the 
other hand, our measurements suggest that the wavevec- 
tor dependence of the electron-phonon coupling consti- 
tutes a general property of covalent crystals. This aspect 
of the coupling can be understood by a simple tight bind- 
ing model on a linear chain: if the lattice displacement 
modifies the overlap between atomic orbitals the bind- 
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FIG. 4: Peak position of the surface resonance a for parallel 
wave vector k\\ = 0.12 A -1 as a function of pump probe delay 
(red marks, left axis), and electronic temperature (black line, 
right axis). 



ing energy of the Bloch states will change differently for 
different wavevectors. Such k dependence of the electron- 
phonon coupling is essential for several problems of solid 
state physics: it accounts for the presence of Charge Den- 
sity Waves that do not respect the nesting condition [T6] 
and it allows for the correct estimate of the high transi- 
tion temperature in superconducting MgE>2 [17] . 

In the case of bismuth, our data offer an unique op- 
portunity to test the deformation potential extracted by 
DFT. We calculate the band structure of bismuth for 
several values of the displacement u. Figure [3^1 shows 
the dispersion of band b for u between 0.1% and 0.3% of 
the body diagonal length c = 11.86 A. In agrement with 
the measurements, the binding energy of band b strongly 
varies near the T point whereas is less sensitive to the 
phonon coordinate for fey = 0.2 A -1 . On the other hand, 
Figure[3p shows that the wavevector dependence of 
is stronger in the calculations than in the experiment. A 
more accurate simulation of the measured signal would 
require the computation of one-step photoemission spec- 
tra [EE]. 

It is clear in Fig. |3p that the dynamics of electronic 
states does not depend simply on the lattice coordinates. 
Besides coherent phonon oscillations, an additional ef- 
fect shifts peak b towards negative values. This pho- 
toinduced shift occurs also in the surface resonance and 



surface states. As shown in Fig. [4j the dynamics of peak 
a nearly follows the electronic temperature. Under simi- 
lar experimental conditions, the photoinduced increase of 
binding energy as been observed also in the surface states 
of Gadolinium [18]. An opposite behaviour has instead 
been reported in the electronic states of materials with 
photoinduced phase transitions p~4j [15] . 

The detailed analysis of the electronic shift as a func- 
tion of k will be part of fore coming work. Nonetheless, 
we notice that the magnitude of this component depends 
on the electronic wavevector and band index. A wavevec- 
tor dependent stiffening of the electronic states in the 
bonding band can be expected because of the reduced ion 
core screening at elevated electronic temperature. How- 
ever, the finite temperature DFT calculations indicate 
that the charge transfer of figure [3j3 accounts at maxi- 
mum for 2 meV to the shift observed in band b. This 
discrepancy may be due to the fact that our bulk model 
of the electronic states does not take into account charge 
redistribution and carrier transport at the surface of the 
sample. Alternatively, the dynamical aspect of the elec- 
tronic screening may require simulations that go behind 
the DFT method Qjg. 

In conclusion, the electronic states of photoexcited bis- 
muth oscillate with the frequency of the A\ g mode. By 
analyzing the photoelectron spectra collected at differ- 
ent angles, we reconstruct the deformation potential of 
specific Bloch states. A quantitative agreement with 
DFT calculations confirms that the photoinduced charge 
transfer accounts for the amplitude of the observed os- 
cillations. At last, we mention that a purely electronic 
effect adds to the temporal evolution of the electronic 
states with the dynamics of the electronic temperature. 
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